Effect of the ly-term for a t — U — W Hubbard ladder 
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Antiferromagnetic and dj.2_y2-pairing correlations appear delicately balanced in the 2D Hubbard 
model. Whether doping can tip the balance to pairing is unclear and models with additional 
interaction terms have been studied. In one of these, the square of a local hopping kinetic energy 
Hw was found to favor pairing. However, such a term can be separated into a number of simpler 
processes and one would like to know which of these terms are responsible for enhancing the pairing. 
Here we analyze these processes for a 2-leg Hubbard ladder. 



The interplay of antiferromagnetism and (1^2 _y2 super- 
conductivity in the 2D Hubbard model remains an open 
question. ^ Weak coupling calculations originally sug- 
gested that doping could drive the ground state from 
an antiferromagnet to a dx^-y-i superconductor. |^ How- 
ever, numerical Monte Carlo calculations have found only 
short range dx2-y2 pairing correlations. |^-|^ This may 
be due to the finite lattice sizes that have been studied, 
the difficulty in attaining low temperature results or pos- 
sibly that the t — U Hubbard model lies just outside the 
superconducting parameter regime. 

One approach to this problem is then to add various 
terms to the basic Hubbard model and see what it takes 
to drive it into a superconduting state. In this spirit, a 
recent Monte Carlo study ||^ added a term Hw , involving 
the square of the local hopping kinetic energy around a 
site, 
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with Ki equal to the local kinetic energy involving site i 
and its near neighbors at i -\- 6, 
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With Hw added to the 2D t-U Hubbard model, the half- 
filled system exhibited a transition from an antiferromag- 
netic phase to a d2.2_j^2-pairing phase at a critical value 
of W . Separating Hw into various pieces, it was found 
that it contained one-electron hopping terms, exchange 
interactions and triplet and singlet four particle scatter- 
ing terms. One would like to understand which of these 
terms or what combination of the terms are responsible 
for enhancing superconductivity. Unfortunately because 
of the fermion sign problem it has not been possible to 
carry out a Monte Carlo calculations for the individual 
terms. However, density matrix renormalization group 
(DMRG) techniques § can be used to study the individ- 
ual pieces of the Hw interaction. Here we describe the 



results of such a study for a 2-leg ladder. For such a sys- 
tem, we can determine the effect of the individual terms 
for both the half-filled and the doped system. As we will 
discuss in the conclusion, it is important to note that the 
half-filled 2-leg ladder has a spin gap which distinguishes 
it from the 2D half-filled Hubbard model. Nevertheless, 
it is instructive to see what effect the various parts of W 
have on the pairing correlations for a ladder. 
We begin with the usual Hubbard Hamiltonian 
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with a one electron hopping kinetic energy and an onsite 
Coulomb interaction U . The sum < ij > is over all pairs 
of nearest neighbors. We will measure all energies in units 
of t. We then add the interaction (|l|) with W positive. 
Monte Carlo calculations for a 2D half-filled system with 
the Hamiltonian 



H^Hr 



H 



w 



(4) 



find a quantum phase transition between an antiferro- 
magnetic Mott insulator and a c?2,2_j,2-wave supercon- 
ducting phase when W is increased to a value of order 
0.35. Q However, the 2D Hubbard model at half-filling 
has an antiferromagnetic ground state while a 2-leg lad- 
der is characterized by a spin gap. js) Thus, as we will 
see, the behavior of a two-leg ladder as W is turned on, 
can be different. 

It is convenient to decompose the interaction Hw as 
follows il 
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with 

Hw, = -AWlY,in^'\ + nil) 
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Hw^ = -W2 cI+Sm'^-^+S' ..'y (6b) 
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Hw,=-W,J2^ls^^s (6d) 
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/\/2 are triplet pair creation op- 
erators, and A]^ = [4,A+s,i - 4,A+sa) /V^isa sin- 
glet pair creation operator. If one sets all the Wi equal 
to W, the original Hyy interaction (|l|) is recovered. Here 
we will examine the effect of the individual terms. 
renormalizes the chemical potential and contains 
next-nearest and next-next-nearest neighbor one-electron 
hopping terms. scatters an onsite singlet to neigh- 

bors sites while Hw^ , which comes with a positive sign, is 
a triplet scattering term. Finally Hw^ and Hwe involve 
singlet pairs. It had been thought for the 2D system that 
the relevant terms for the quantum transition were -ffivs 
and Hwe ■ i] 

Here, in order to determine the effects of the individual 
terms, we have studied the model on a two-leg ladder us- 
ing DMRG techniques. All the runs were done on 2 x 32 
ladders keeping up to 800 states leading to a maximum 
discarded weig ht of 10"*^. We calculated the singlet pair- 
ing correlation function Dafj{i) defined as 



D,y{£) = (A,(z + £)At(z)) 

Dyy{£) ^ {Ayit + i)Al{z)) 



(7) 

(8) 
(9) 



where Aa{i) 



^1,14+5, 



- 4,i4+s„.p = (I'O) and 
dy = (0, 1). For clarity, in the following we show the rung- 
rung correlation function Dyy{£). D^xi^) and Dyy{£) 
were always positive while Dxy{t} was negative corre- 
sponding to a d^2_y2-\ik& strucure. 

The results for the half-filled case with J7 = 4 and 
Wi = or 0.25 are shown in Fig. |l|. In the plot of 
Dyy{£) we have kept £ < 12, with the measurements made 
in the central portion of the ladder. In this region the 
effects of the open ends are negligible. We clearly see 
in part (a) that when all Wi are turned on there is an 
enhancement of the pairing (as found in the 2D Monte- 
Carlo simulations). However, if we only turn on W^, 
there is a suppression of pairing. For the 2-leg ladder, this 
can be understood by noting that Hwr, can be written as 
an antiferromagntic exchange interaction 



Hw, = 2W; 



(10) 



Now as one knows, ||] a 2-leg Heisenberg ladder has a 
spin gap As « 0.51J. Thus the effect of at half- 
filling is to increase the spin gap by a factor of order W5 
and this leads to an exponentially more rapid decay of 
the pairing correlations. 
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FIG. 1. The rung-rung singlet pairing correlation function 
Dyy{£) versus ^ on a half- filled ladder with U = A for (a) all 
W^ = 0, only Wi, = 0.25 and all W« = 0.25 and (b) all = 
and only one Wi — 0.25. 

Fig. (b) shows the effect of the other terms. Hwi 
has no effect as expected since it just renormalizes the 
chemical potential and we have fixed (n) = 1. The ad- 
ditional one electron hopping term leads to only a 
small change in the pairing. However which scatters 
an onsite singlet to neighbors sites enhances the pairing 
despite the presence of U which lowers the double occu- 
pancy. 
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FIG. 2. The rung-rung singlet pairing correlation function 
Dyy{l) versus ^ on a half-filled ladder with W3 — 0.25 and 
various values of U. 

We have performed other calculations including only 
which show that for larger U this enhancement is 
supressed as one would expect (see Fig. ^). Neverthe- 
less, for U/t = 4 where the previous 2D Monte Carlo 
calculations were run, contributes to enhanced the 
pairing. iJvKi also leads to enhanced singlet pairing. Note 
that it has a positive coefficient which suppresses triplet 
pairing, leaving more phase space for singlet pairing. Fi- 
nally Hwei which describes singlet pair hopping for (i,5) 
to {i,S'), also enhances the pairing. We should point 
out that although there is an enhancement in the pair- 
ing, this enhancement is in fact relatively small for the 
two-leg ladder for reasonable values of W and as U is in- 
creased, so that the W3 term is reduced, the suppression 
of the pairing by the W5 term becomes dominant. This is 
clearly seen in Fig. ^, where we show the pairing correla- 
tion function for the half-filled ladder with all Wi = 0.25 
and various values of U . 
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We now investigate the nature of the magnetic ordering 
by examining the spin-spin correlation function 

S{r) = (S+Si^^) (11) 

where 5'/ {S^) are the spin raising (lowering) operators 
corresponding to — c\^ass'Ces- We then perform a 
Fourier transform to obtain the static structure factor 

5(g) = ^e'^'-5(r). (12) 

r 

The resulting structure factor is plotted in Fig. ^ for var- 
ious Wi interactions. As previously noted, the effect of 
W5 is to increase the spin gap which leads to a broaden- 
ing of the Lorentzian and a decrease of its peak at (tt, tt). 
When all of the Wi terms are present, S{q) appears to 
simply be reduced in magnitude at all q values indicat- 
ing a reduction of the local moment \/JSf) due to the 
delocalization effects of Hw- 
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FIG. 4. Structure factor S{q) of the spin-spin correlation 
function versus for the half-filled system and qy — tt. The 
inset shows an enlargement of the region near {n,n). 

We now turn to the doped case and consider the same 
lattice with [/ = 4 and 8 holes corresponding to a filling 
(n) = 0.875. Fig. | show results for Dyy{t). We clearly 
see in part (a) that in this case the inclusion of -ffv^s 
enhances the pairing while in part (b) we see that all of 
the remaining terms are essentially irrelevant. Thus the 
14^5 term, which corresponds to adding a near neighbor 
exchange J = 2W5 enhances the pairing correlation in 
the doped system. 



FIG. 3. The rung-rung singlet pairing correlation function 
Dyy{£) versus ^ on a half-filled ladder with all Wi = 0.25 and 
various values of U. 
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FIG. 5. The rung-rung singlet pairing correlation function 
Dyy{£) versus ^ on a doped ladder with (n) — 0.875 and U — 4 
for (a) all Wi = 0, only W5 = 0.25 and aU Wi = 0.25 and (b) 
all Wi=0 and only one Wi = 0.25. 

Thus we conclude, that while Hw with [/ = 4t can 
slightly enhance the pairing correlations of a half-filled 
ladder, this is in fact a small effect. Furthermore, for 
large values of U/t, Hw leads to a suppression of the 
half- filled pairing correlations. This can be understood 
in terms of the dominance of ffvi/j , which represents an 
effective antiferromagnetic exchange increasing the spin 
gap and suppressing the pairing correlations. However, 
for the doped ladder, W5 acts to enhance the pairing 
correlation since it increases the effective exchange inter- 
action and the pair binding energy. Clearly, in light of 
the 2D results, where it was found that Hw could lead 
at half-filling to a dx2-y2 pairing state, one would like to 
extend the DMRG calculations to a 3-leg ladder which 
has a vanishing spin gap at half-filling. 
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